EPJ manuscript No. 

(will be inserted by the editor) 



o 
o 

(N 



Ph. 



Analysis of X-ray spectra emitted from laser-produced plasmas 
of uranium 

J. P. Marques 1 , P. Indelicato 2 , C. I. Szabo 2 , and F. Parente 3 

1 Centro de Fi'sica Atomica e Departamento Fi'sica, Faculdade de Ciencias, Universidade de Lisboa, 
Campo Grande, Ed. C8, 1749-016 Lisboa, Portugal, e-mail: jmmarques@fc.ul.pt 

2 Laboratoire Kastler Brossel, Ecole Normale Superieure; CNRS; Universite P. et M. Curie - Paris 6 

Case 74; 4, place Jussieu, 75252 Paris CEDEX 05, France, e-mail: paul.indelicato@spectro.jussieu.fr 

3 Centro de Fisica Atomica da Universidade de Lisboa e Departamento Fi'sica da Faculdade de Ciencias e Tecnologia da 
Universidade Nova de Lisboa, Monte da Caparica, 2825-114 Caparica, Portugal, e-mail: facp@fct.unl.pt 

Received: January 19, 2009/ Revised version: date 



o 

H— > 

cd 

c/a 
O 



> 
On 
t> 
l> 
(N 

o 

OS 

o 



Abstract. In this paper, we used the multiconfiguration Dirac-Fock method to generate theoretical X-ray 
spectra for Co-, Ni-, Cu-, Zn-, Ga-, Ge-, As-, Se-, Br-, Kr-, and Rb-like uranium ions. Using the distribution 
of these ions in a laser-produced plasma, for different plasma temperatures, we generate theoretical spectra 
which are compared to experimental data. 

PACS. 32.30.Rj - 31.15.A - 32.70.Cs 



1 Introduction 

Hard X-ray spectra from laser-produced plasmas have 
been the subject of extensive research in the last few years. 
This is mainly due to their importance in what concerns 
alternative energy sources and atomic, plasma and mate- 
rials physics research, including astrophysics. In a recent 
review Hudson et al. lj emphasize the importance of this 
subject for laser-driven inertial confinement fusion, exten- 
sively describing their challenges, the efforts that were un- 
der way in the world, and suggesting possible responses. 

Usually, in this kind of plasmas a large number of 
atomic species are present, in different stages of ionization, 
and, within each of these stages, in different energy lev- 
els, the number of the latter being sometimes very large. 
Emitted X-radiation may be an important diagnostic of 
the plasma provided that one is able to connect the de- 
tected radiative lines with the emitting ion. Therefore, a 
calculation as precise as possible of the energy levels of the 
ions which may be present in a plasma is of crucial im- 
portance for the interpretation of recent high-resolution 
experimental spectra. 

Seely et al. [2] reported on X-ray spectra from laser 
produced plasmas obtained when planar foils of U and Pb 
were bombarded with a total energy of about 12 kJ at the 
OMEGA (University of Rochester's Laboratory for Laser 
Energetics, LLE) laser facility. The spectra show a few 
intense and narrow features in the 12-22 keV energy range, 
interpreted as inner-shell transitions resulting from L-shell 
vacancies created by energetic electrons. The calculations 
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of these authors suggest that the transitions are Ni-like or 
lower ionization stages. 

Yuan and Moses [3] developed an atomic code, YAC, 
to calculate non-LTE plasma spectra using the detailed 
accounting method for line transitions. With this code 
they obtained the ionic distributions for different plasma 
temperatures from 450 eV to 850 eV. Comparing their re- 
sults with experimental Seely et al. [2] uranium spectra, 
they estimate that plasma temperature was 650 eV. The 
plasma conditions were estimated with the help of the ID 
hydrodynamic code BUCKY g). 

In this work we use the Multi-Configuration Dirac- 
Fock code of Desclaux and Indelicato [5j[6] , which includes 
QED corrections, to compute the energies of several X- 
ray lines from Co-like to Rb-like uranium ions which are 
important for the interpretation of Seely et al. uranium 
spectra. 

Our results, together with the ionic distributions of 
Yuan and Moses allow for a more precise determination 
of the plasma temperatures in Seely et al. [2] experiment. 



2 Relativistic calculations 

We calculated bound-states wave functions using the 2007 
version of the Dirac-Fock program of J. P. Desclaux and 
P. Indelicato [6] . Details on the Hamiltonian and the pro- 
cesses used to build the wave-functions can be found else- 
where [5117115119]. 

The total wave function is calculated with the help of 
the variational principle. The total energy of the atomic 
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system is the eigenvalue of the equation 
H na JjM (. . . 7 nj . . .) = Enj.M&nj.Aii- ■ ■ ,n, . . .), 

(1) 

where 77 is the parity, J is the total angular momentum 
eigenvalue, and M is the eigenvalue of its projection on 
the z axis J z . In this equation, the hamiltonian is given 

by 

JV 

H"° pai ' = £ Hd W + E V U ( I r « I ) ' ( 2 ) 

i— 1 i<j 

where He is the one electron Dirac operator and V« is an 
operator representing the electron-electron of order one in 
a. The expression of Vij in Coulomb gauge, and in atomic 
units, is 



Vij = 



1 

OL t ■ OLj 



(3a) 
(3b) 



^[cos(^) 1 



cos (^iilil) - i 

+ c 2 (a t ■ V<)(oy • V,) , (3c) 

u} ij r ij 

where = — rj \ is the inter-electronic distance, Uij is 
the energy of the exchanged photon between the two elec- 
trons, a.i arc the Dirac matrices and c is the speed of light. 
We use the Coulomb gauge as it has been demonstrated 
that it provides energies free from spurious contributions 
at the ladder approximation level and must be used in 
many-body atomic structure calculations [lOilllj . 

The term (|3a| represents the Coulomb interaction, the 
term (|3b[) is the Gaunt (magnetic) interaction, and the last 
two terms (|3c|) stand for the retardation operator. In this 
expression the V operators act only on and not on the 
following wave functions. 

By a series expansion of the operators in expressions 
(|3bp and (|3c[) in powers of ujijrij/c <C 1 one obtains the 
Brcit interaction, which includes the leading retardation 
contribution of order 1/c 2 . The Breit interaction is, then, 
the sum of the Gaunt interaction (|3b[) and the Breit re- 
tardation 



B 



Oil ■ OLj {OLi ■ Vij) (dij ■ Vij) 



2r„, 



2r 3 - 



(4) 



In the many-body part of the calculation the electron- 
electron interaction is described by the sum of the Cou- 
lomb and the Breit interactions. Higher orders in 1/c, de- 
riving from the difference between Eqs. (f3"c)) and ([U are 
treated here only as a first order perturbation. All calcula- 
tions are done for finite nuclei using a Fermi distribution 
with a thickness parameter of 2.3 fm. The nuclear radii 
are taken from reference [12J . 

The MCDF method is defined by the particular choice 
of a trial function to solve equation ([I]) as a linear combi- 
nation of configuration state functions (CSF): 



n,j,M) = £ Cv \ y > n > J > M ) ■ 



(5) 



The CSF are also eigenfunctions of the parity 77, the total 
angular momentum J 2 and its projection J z . The label v 
stands for all other numbers (principal quantum number, 
...) necessary to define unambiguously the CSF. The c„ 
are called the mixing coefficients and are obtained by di- 
agonalization of the Hamiltonian matrix coming from the 
minimization of the energy in equation (|5|) with respect 
to the c„. 

The CSF are antisymmetric products of one-electron 
wave functions expressed as linear combination of Slater 
determinants of Dirac 4-spinors 



\v,n,j,M) = ^di 



(C) 



where the ip-s are the one-electron wave functions and the 
coefficients di are determined by requiring that the CSF is 
an eigenstate of J 2 and J z . The di coefficients are obtained 
by requiring that the CSF are eigenstates of J 2 and J z . 
The one-electron wavefunctions are defined as 



tp(r) 



ij£*(fl)Q(r) 



(7) 



where xt i s a two-component spinor, and P and Q are 
respectively the large and small components of the wave- 
function. A variational principle provides the integro-dif- 
ferential equations to determine the radial wave functions 
and a Hamiltonian matrix that provides the mixing coef- 
ficients c„ by diagonalization. One-electron radiative cor- 
rections (self-energy and vacuum polarization) are added 
afterwards. All the energies are calculated using the ex- 
perimental nuclear charge distribution for the nucleus. 

The so-called Optimized Levels (OL) method was used 
to determine the wave function and energy for each state 
involved. This method allows for a full relaxation of both 
initial and final states providing much better energies and 
wavefunctions. However, spin-orbitals in the initial and 
final states are not orthogonal, since they have been op- 
timized separately. The formalism to take in account the 
wave functions non-orthogonality in the transition proba- 
bilities calculation has been described by Lowdin [13j . The 
matrix element of a one-electron operator O between two 
determinants belonging to the initial and final states can 
be written 

(v>nJM\ E O {n) \v'n'J'M') = 



Nl 



ipi (r m ) ■■■ip m (r m ) 

0i (ri) ■■■ 4> m (n) 



4>i (r m ) ■ ■ ■ 4> m (r m ) 



(8) 



where the ipi belong to the initial state and the 4n and 
primes belong to the final state. If ip = \nnfi) and 4> — 
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|nV//) are orthogonal, i.e., (n,K[i\n' k' //) = 5 
the matrix element ([5]) reduces to one term (ipi\0\(j>i) 
where i represents the only electron that does not have 
the same spin-orbital in the initial and final determinants. 
Since O is a one-electron operator, only one spin-orbital 
can change, otherwise the matrix element is zero. In con- 
trast, when the orthogonality between initial and final 
states is not enforced, one gets [T3] 



(viijm\ o (n) yn'j'W) = 

i=l 



(9) 

where D^i is the minor determinant obtained by crossing 
out the ith row and j'th column from the determinant of 
dimension N x N, made of all possible overlaps (ipk\<l>i')- 
Radiative corrections are also introduced, from a full 
QED treatment. The one-electron self-energy is evaluated 
using the one-electron values of Mohr and coworkers [T31 
[TMTB] and corrected for finite nuclear size [T7]. The self- 
energy screening and vacuum polarization are treated with 
an approximate method developed by Indelicato and co- 
workers plfTOpOpT]- 



3 Results and discussion 

In the analysis of their laser-produced uranium spectrum, 
using the HULLAC code [22], Seely et al. 2 concluded 
that the observed energies of the two most intense fea- 
tures in the spectrum are 260 and 230 eV higher than the 
neutral U Lai and La2 energies, respectively. Other fea- 
tures were observed at energies 90 and 260 eV higher than 
the L/3i energy. These energy differences led the authors 
to the conclusion that the observed transitions were due 
to highly charged ions in the vicinity of the Ni-like U or 
in lower ionization stages, around As- and Si-like U. 

Yuan and Moses [3] analysis of the same spectrum by 
means of an atomic model (YAC) which uses the detailed 
configurations as effective levels for the population equa- 
tions but uses the detailed term accounting method for the 
lines transitions. Relativistic formulas are adopted in order 
to apply to high-Z elements. A comparison of the oscilat- 
tor strengths and energies between YAC and GRASP 23J 
results for Cu-like U 2p-4d transition shows that although 
the oscillator strengths agree very well, differences of 40 
to 80 eV in the transition energies were found. Further- 
more these authors present ionic distributions for different 
plasma temperatures from 450 to 850 eV and estimate, 
from their analysis, that the plasma is of a temperature 
of 650 eV. 

In this work we calculated the energies of all levels in 
the U (t = 1,2,3) Mi (i = 1, • • • , 5), N, (i = 1, • • • , 7) 
and Oi (for Rb only) one hole configurations in Co, - 
Ni-, Cu-, Zn-, Ga-, Ge-, As-, Se-, Br-, Kr- and Rb-like 
uranium. All transition energies and probabilities between 
initial Lj and final Mj, Nj and Oi states were computed 
taking into account the interaction of the inner holes with 
electrons in outer unfilled shells. In this calculation no 
electronic correlation was taken into account, but fully 



relaxed initial and final wavefunctions were independently 
calculated, so full relaxation in both energy and transition 
rates is included. 

The intensity of the line corresponding to a transition 
from level i to a level j in a U 9+ ion with and in level i, 
is given by 

1% = hvijAfjN?. (10) 

In this equation hvij and A^- are the radiative i —* j 
transition energy and probability, respectively, and Nf is 
the density of U 9+ ions in level i. This density is obtained 
from the balance equations [53] and equation (|10[) can be 
written as 



hv. : 3 



N? 



A '■' 



-U) 



(11) 



Here A\ is the transition for deexcitation by all possible 
processes (radiative and radiationless), u>i is the fluores- 
cence yield for level i and the sum Y]^, N? is performed for 
all charged states q' . In this work the fluorescence yields 
given by Krause [35] for the neutral atoms were used. 

In the case of Zn-like (ls 2 2s 2 2p 6 3s 2 3p 6 3d 10 4s 2 ) and 
Kr-like (ls 2 2s 2 2p 6 3s 2 3p 6 3d 10 4s 2 4j? 6 ) ions we calculated 
the energies and transition probabilities for all lines. The 
results are given in table [T] 

For the other charge states we have to account for the 
interaction of the inner hole with the electrons in the un- 
filled outer shells. For example, for the Is 2 2s 1 2p 6 3s 2 3p 6 
3d 8 4s 2 configuration in Ni-like uranium with one hole in 
Li subshell, J = 1/2, • • • ,9/2, leading to 20 energy lev- 
els. In the case of an hole in L2 or L3 subshells, i.e., the 
Is 2 2s 2 2p 5 3s 2 3p 6 3d 8 4s 2 configuration, J values range 
from 1/2 to 11/2 leading to 45 levels. So, for the initial 
states we have a large number of levels, 65 in this case. 
In the final Mi to Ni states a similar situation will oc- 
cur. The obvious consequence is that for a full calculation 
we have to compute a large number of transitions. In this 
case we calculated all transition probabilities between ini- 
tial and final levels for Aj = 0, ±1 which represents 3234 
lines. A similar situation occurs for all the other charge 
states that were calculated, except Zn-like and Kr-like ura- 
nium. The transition probabilities and energies calculated 
in this work can be found in [26] . 

In Fig. [T] the theoretical spectrum obtained using for 
the lines a Gaussian distribution with a width of 95 eV, 
is presented for Ni-like uranium. In the inset one can see 
a detailed view of the Lai^ (L3— M^s) structure. 

In Fig. [2] we compare the calculated spectra for Co-, 
Ni-, Cu-, Zn-, Ga-, Ge-, As-, Se-, Br-, Kr-, and Rb-like 
uranium. In this figure it is clear that as Z increases, a 
shift of the peak energies to the lower energy side of the 
spectrum is observed. 

In order to compare with the experimental spectrum 
obtained by Seely et al. [2] we present in Fig. [3] the theo- 
retical spectra calculated for plasma temperatures of 650, 
750, and 850 eV, for a Gaussian width of 95 eV, plotted 
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Fig. 1. Calculated Ni-like uranium X-ray spectrum. This figure was generated using a Gaussian width of 95 eV. In the inset is 
presented a detailed view of the Lai, 2 structure. 
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Fig. 2. Calculated spectra for Co-, Ni-, Cu-, Zn-, Ga-, Ge-, As-, Se-, Br-, Kr-, and Rb-like uranium. 
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Table 1. X-ray lines calculated energies (in eV) for Zn- and Kr-like uranium. 





initial 




final 


Zn-like 


Prob. 


Kr-like 


Prob. 


LI 


2Si/2 


Ml 


3Sl /2 


16430.795 


8.49E+10 


16348.010 


8.29E+10 


LI 


2si /o 


M2 


3r>i /i 
°Fl/2 


16782.763 


1.38E+15 


16702.777 


1.35E+15 


LI 


2si /'-> 

io l / 2 


M3 


3Dq /o 


17708.851 


1.25E+15 


17616.946 


1.21E+15 


LI 


2si 

iJ l / 2 


M4 


3d3/2 


18234.176 


6.62E+13 


18162.897 


6.49E+13 


LI 


2Si/2 


M5 


3ds/2 


18422.600 


9.94E+13 


18347.679 


9.73E+13 


LI 


2Si/2 


Nl 


4Si/2 


21095.510 


5.31E+10 


20866.003 


4.95E+10 


LI 


2si 


N2 


4r>i 

^^1/2 






21036.126 


4.68E+14 


LI 


2si /'-> 


N3 


4Dq /o 






21340.323 


5.14E+14 


L2 


2Di /o 


Ml 


3si /o 


15631.985 


1.57E+14 


15539.004 


1.54E+14 


L2 


2Di /o 


M2 


3di /o 

U F1 / 2 


15983.953 


2.19E+10 


15893.772 


2.12E+10 


L2 


-^Pl/2 


M3 


3r>Q 

°Fd/2 


16910.041 


8.31E+12 


16807.940 


7.94E+12 


L2 


2di /o 


M4 


3d-s /o 


17435.366 


5.47E+15 


17353.891 


5.36E+15 


L2 


2pi/2 


M5 


3d5/2 


17623.790 


6.95E+11 


17538.673 


6.73E+11 


L2 


2Pl/2 


Nl 


4Sl/2 


20296.700 


5.29E+13 


20056.998 


5.12E+13 


L2 


2P1/2 


N2 


4pi/ 2 






20227.120 


1.20E+10 


L2 


2P1/2 


N3 


4p 3 /2 






20531.317 


3.27E+12 


L3 


2P3/2 


Ml 


3Sl/ 2 


11829.419 


2.76E+14 


11749.368 


2.71E+14 


L3 


2P3/2 


M2 


3P1/2 


12181.388 


3.59E+12 


12104.136 


3.48E+12 


L3 


2P3/2 


M3 


3P3/2 


13107.475 


3.39E+12 


13018.305 


3.26E+12 


L3 


2P3/2 


M4 


3d 3/2 


13632.801 


4.59E+14 


13564.256 


4.50E+14 


L3 


2P3/2 


M5 


3d 5/2 


13821.225 


4.06E+15 


13749.038 


3.98E+15 


L3 


2P3/2 


Nl 


4si/ 2 


16494.135 


8.75E+13 


16267.362 


8.41E+13 


L3 


2P3/2 


N2 


4pi/ 2 






16437.485 


9.48E+11 


L3 


2P3/2 


N3 


4P3/2 






16741.682 


1.17E+12 



against the experimental data. For this purpose, we used 
the ionic distributions for the different plasma tempera- 
tures, taken from Yuan and Moses [3]. 

We conclude that the calculation for 850 eV plasma 
temperature is the one that best fits the experimental 
data, or, comparing the plots for 850 eV and 750 eV we 
may guess that the plasma temperature should be between 
these two values. Trying to verify this conclusion, we per- 
formed a x 2 test to see if the spectrum obtained with one 
of the temperatures matches better than the others the ex- 
perimental data. The best matches were obtained for 750 
and 850 eV which are close to each other, nevertheless 
the best match was obtained for a plasma temperature of 
850 eV. This is in contrast with Yuan and Moses [3] con- 
clusions that based in their own calculations pointed to a 
plasma temperature of 650 eV. 

As in the calculations by other authors the major dis- 
agreement with the experiment can be seen in the intensity 
ratio between the Levi and L«2 lines. This behavior has 
been explained by opacity and plasma density gradient ef- 
fects which may have an important role on the shape of 
the observed spectrum [2]. 
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